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Our stress relaxation measurements on wormlike micelles using a Rheo-SALS (rheology + small 
angle light scattering) apparatus allow simultaneous measurements of the stress and the scattered 
depolarised intensity. The latter is sensitive to orientational ordering of the micelles. To determine 
the presence of causal influences between the stress and the depolarised intensity time series, we 
have used the technique of linear and nonlinear Granger causality. We find there exists a feedback 
mechanism between the two time series and that the orientational order has a stronger causal effect 
on the stress than vice versa. We have also studied the phase space dynamics of the stress and the 
, depolarised intensity time series using the recently developed technique of cross recurrence plots 

' (CRPs). The presence of diagonal line structures in the CRPs unambiguously proves that the two 

time series share similar phase space dynamics. 
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Wormlike micelles under shear flow are turning out to be classic examples of chaotic nonlinear dynamical systems. 
Wormlikc micelles are long, flexible cylindrical objects formed by the self-assembly of surfactant monomers. The 
' individual worms entangle above a critical concentration to form viscoelastic gels. Wormlike micelles relax stress 
through two mechanisms: reptation, analogous to polymers, and by scission and recombination Due to the 
latter process of relaxing stress, wormlike micelles show extreme shear-thinning behaviour and this shows up as a flat 
plateau in the flow curve (stress vs. shear rate). Gates et al. were the flrst to predict this kind of flow behaviour 
and they attribute this to a mechanical instability of the shear banding type In shear banding systems, the 
system splits into coexisting bands that support the same stress but have different average shear rates. The high 
shear rate band is lower in viscosity and is the nematic phase, while the low shear band has a higher viscosity and 
■ is the isotropic phase. Stress relaxation measurements in the plateau region of the flow curve for viscoelastic gels 
of surfactant cetyltimethyammonium tosylate (GTAT) showed chaotic dynamics ,3:|. Inspite of the short data trains 
OO ' analysed (w 5000 data points), owing to the low dimensional nature of the attractors, the phase space dynamics 
could be quantifled in terms of Lyapunov exponents and fractal correlation dimensions. This phenomenon, termed 
. rheochaos, has now been observed in many other systems d, d, 0, 0] and recently the route to rheochaos has also 
y—i ' been observed. Ganapathy and Sood have shown that the route to rheochaos is via Type-II intermittency in stress 
C , relaxation measurements Q for wormlike micellar gels of surfactant GTAT in the presence of salt sodium chloride 
(NaGl). Depolarised small angle light scattering measurements done simultaneously with flow experiments showed 
^ . that the depolarised scattered intensity, sensitive to orientational order fluctuations, showed similar dynamics as that 
JJ"! ■ of the stress/shear rate [l). 

rS ' In the experiments reported so far 0, 0, spatial inhomogeneity always accompanies rheochaos. Spatial inhomogene- 
^ . ity has been captured in many recent theoretical models of rheochaos [l^, El, ■ The theoretical model relevant 

" " ' to the work described in this paper, proposed by Ghakrabarti et al. |12l |. is a spatial extension of the homogeneous 
model first proposed by Hess et al. [13]. This model treats the system as an orientable fluid and considers the 
spatio-temporal evolution of the symmetric traceless nematic order parameter . The equations of the motion of 
the stress are not fundamental but are derived from the underlying dynamics of the alignment tensor or the local 
orientational order parameter. An important assumption here is that the orientational order parameter equations are 
studied in the passive advection approximation i.e. the effect of the stress field, arising from orientational ordering, 
on the fiow profile is ignored. Although this model ignores features specific to wormlike micelles like breaking and 
recombination and shear induced elongation or breaking of the micelles, it captures the essential orientable nature of 
the worms under flow. This model also predicts the route to rheochaos to be via spatiotemporal intermittency (l^ . 

Given this strong theoretical backing, it would indeed be useful to check if there exists any cause and effect 
relationship between the time series of stress and the depolarised scattered intensity that were measured experimentally 
for the GTAT + NaGl system Q . If the passive advection approximation used in the theoretical model [12] described 
above is indeed justifled, one should expect to find that the orientational order has a stronger causal influence on 
the stress than vice versa. We use the technique first developed by Granger to study linear causal effects [3] and 
a recent extension of it for nonlinear time series [l5| to determine if there exists any causality effects in our data. 
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In this paper we show that there exists a Granger feedback mechanism between the two time series and that the 
depolarised intensity or nematic ordering has a stronger influence on the stress than vice versa. Secondly, since 
standard nonlinear time series analysis procedures fail to quantify the dynamics of the chaotic region in intermittent 
data due to insufficient data points, we have used the method of cross recurrence plots (CRPs). We have analysed a 
data train that shows Typc-II intermittency and wc find that the stress and orientational order (depolarised intensity) 
share similar phase space behaviour in the laminar and chaotic regions. 



II. THEORY 



A. Linear Granger Causality 

Consider an experiment in which two time series X{t) and Y{t) have been simultaneously measured. If the prediction 
of X{t) at the current time can be improved by incorporating past information from Y{t), then Y(t) causes X{t). 
This idea of causality was formalised by Granger with reference to linear regression models of stochastic processes 
(T^ and has numerous applications in economics [3 ■ The study of causality has been also applied to medicine [TtI 
and nonlinear dynamical systems p^ . Quantitatively speaking, Y{t) has a causal influence on X{t), if the variance 
of the prediction errors of X{t) decreases on including information from Y(t). The basis for Granger causality lies in 
linear prediction theory. Let X{t) be a stationary process. It can be expressed using an autoregressive (AR) model 
where the current value of X{t) depends on m past values. 

m 

X{t)^J2a,X{t~i) + ex{t) (1) 

i=l 

Here ex(i) is the error in prediction whose magnitude can be evaluated by its variance var(ex(0) ^^'^ o^j's are the 
regression coefficients. If we now include information from another time series, Y{t), that has also been simultaneously 
measured, then the prediction of the current value of X{t) based on its own past and the past values of Y{t) is 

m m 

X{t) = «»^(^ ^'^(^ - *) + ^x\Y{t) (2) 

1=1 i=l 

Here ex|y(i) is the prediction error in X{t) given information from Y{t) and a^, hi are the regression coefficients. If 
var(ex|y(i)) < var(ex(^)), then Y{t) has a causal influence on X{t). Similarly, for Y{t) 

m 



m m 

Y{t) = - + I] d,X{t eY\x{t) (4) 

i=l i=l 

If var(ey|x(i)) < var(ey(t)), X{t) has a causal influence on Y{t). If var(ey|x(i)) < var(ey(t)) and var(ex|y(0) < 
va.r{ex{t)), then there exists a feedback mechanism between X{t) and Y{t). 



B. Extended Granger Causality 

Granger causality was developed for linear processes and a direct application on nonlinear systems might not be 
appropriate. Various theoretical models have been put forward to extend Gr ang er causality for nonlinear systems. 
We briefly describe below the steps involved in a recent model by Chen et al. [15j . 

Consider two time series X{t) and Y{t). Their joint dynamics can be reconstructed using the following delay vector 

MM 

z(t) = (x(tr,y(tn (5) 



where 

x(t) = (X(i),X(t-ri),...,X(t-(mi-l)Ti))^, 
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y{t) = {Y{t),Y{t — Ti),...,Y{t — (m2 — 1)t2))'^, mi, 1712 are the embedding dimensions for the two time series 
determined from the method of false nearest neighbors, Ti,T2 are the delays determined from the first minimum of 
the mutual information function and the superscript T stands for transpose. To investigate causal influences the flow 
of time is important and we set ti = T2 = t. 

In the reconstructed joint phase space of z(t), consider a small local neighborhood of size e. The dynamics of 
this local neighborhood can be described within a linear approximation and a linear autoregressive model can be 
used to predict the dynamics within the neighborhood. The errors in prediction are given by tY\x{t) ^-nd tx\Y{i)- 
The reconstruction and the fitting procedure are now performed on the individual time series X{t) and Y{t) in the 
same local neighborhood and the errors ex{t) and tvit) are determined. The error ratios Y&r{eY\x{t)) l^3.v{eY{t)) and 
var(ex|y(i))/var(ex(i)) are computed. Since for a nonlinear system, causal inferences cannot be made by studying 
a single local neighborhood, the entire process described above is repeated for various regions on the attractor and 
the average of the error ratios are computed. The extended Granger causality index (EGCI) is defined as Ay^x = 
(1 — var(ex|r(^))/var(ex(i))), where (•) stands for the average. The EGCI is now computed as a function of the 
neighborhood size e. For linear systems this index will remain roughly constant as e decreases, whereas for nonlinear 
systems, in the small e limit, the EGCI reveals the true nature of the causal infiuence. 

C. Cross Recurrence Plots 

A major drawback of standard nonlinear time series analysis routines is the requirement of a long time series. One 
of the techniques developed to overcome this limitation is the method of Recurrence Plots (RPs) [2l|, [22, [23] . RPs is 
a tool used to visualise phase space trajectories and as the name suggests RPs represent the recurrence of trajectories 
in phase space, a basic feature of nonlinear dynamical systems. The first step towards generating RPs is to completely 
unfold the attractor. Using the Taken's embedding theorem, the phase space trajectory, x(t), can be reconstructed 
by time delayed vectors of the time series X{t), 

x(t) = .t1 = {X{t),X{t - n), ...,X{t ~{m~ 1)t)). (6) 

Here m and t are the embedding dimension as described in the previous section. The recurrence plot is defined as 

R,,, =e(e,-||xl-x;||) (7) 

Here ei is a suitably defined cut-off distance, is the Euclidean norm and 8 is the Heaviside function. If Xj falls 
within distance of Xi, then the two trajectories are close by in phase space and R^j- = 1, else R^j- = 0. The I's and 
O's are color coded and represented as a plot. 

A recent extension of RPs to bivariate time series is the method of Cross Recurrence Plots (CRPs) [13, [1^ [1^ ■ 
Here both the time series are simultaneously delay-embedded in the same phase space and a test for closeness of the 
two trajectories is carried out. The CRP is defined as 

CR,,, =e(e.-||x*-2/;-||) (8) 

RPs yield a wealth of visual information. Trajectories that share similar local dynamics show up as diagonals. 
Static regions in the trajectory are seen as vertical or horizontal structures. In order to be more quantitative, Zilbut 
and Webber have proposed the following recurrence measures [l^, [2^. 1. %Rccurrence: This is defined as the % of 
recurrence points falling within a given radius. 2. %Determinism: This is defined as the number of recurrence points 
forming diagonal line structures. The more %determinism, the more predictable is the signal. Stochastic signals have 
no predictability, whereas chaotic signals have short-term predictability 3. MaxLine: This is the length of the longest 
diagonal. This is inversely related to the maximum Lyapunov exponent which measures the rate at which phase space 
trajectories diverge. 4. Entropy: This is the Shannon Information entropy and is a measure of randomness of the 
signal. It has units of bits/bin. A periodic signal has zero entropy since all the diagonal are of the same length. A 
random process has a finite non-zero Shannon entropy. 5. Trend: This is a measure of the stationarity of the time 
series. 6. %Laminarity: This is a measure of the % of recurrence points forming vertical or horizontal line structures. 
7. Trapping Time: This is a measure of the average length of the vertical or horizontal line structures. The last two 
measures were proposed by Marwan et al. (27j . 

III. EXPERIMENTAL METHODS 

The data analysed in this paper were obtained from our earlier stress relaxation experiments on the wormlike 
micellar system of CTAT 2wt% -I- lOOmM NaCl described in detail in Q. Under controlled shear rate conditions 
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this system shows the Type-II intermittency route to rheochaos. We briefly describe here the method that was 
used in our earher work for generating the bivariate time series of stress and depolarised intensity. The experiments 
were performed on a MCR 300 stress controfled rheometer with SALS attachments at a temperature of 26.5°C. The 
experiments were carried out in a cyhndrical Couette geometry with top and bottom windows made of quartz glass 
(inner cylinder diameter 32mm, height 16.5mm and gap 2mm). A vertically polarised (V) laser beam (A — 658nm and 
spot size 1mm) enters the gap between the cylinders along the vorticity (V x v) direction, where v is the velocity field. 
An analyser below the Couette geometry allows us to select either the vertically (referred as VV) or the horizontally 
polarised/depolarised scattered light (referred as VH) from the sample without disturbing the measurements. A 
condenser beneath the analyser collects the scattered light dominantly from a plane 6mm above the bottom plate and 
forms the image on a screen in the (v, Vv) plane. A 8-bit colour CCD camera (Lumenera 075C, 640 x 480 pixels, 
maximum frame rate - 60fps) at a frame rate of 1 frame/750ms imaged the screen. About 3000 images were grabbed 
for each polarisation while stress relaxation measurements were simultaneously going on. The intensity at various 
wavevectors from the noise filtered image was measured and a time series was generated by repeating the process over 
each image. The data analysed in this paper are at a fixed wave vector, q = 0.75/xm~^. 

IV. RESULTS 

Figure [1^ & b show the partial time series of the stress and the VH intensity at a fixed shear rate of 25 
respectively. The time series shows intermittent behaviour characterised by the presence of laminar and chaotic 
regions. The reason for analysing this particular data set is two-fold: 1. it would be worthwhile to check whether 
the two time series share similar phase space dynamics in the laminar and chaotic regions and 2. the chaotic region 
of the time series cannot be analysed using normal nonlinear time series analysis procedures due to the relatively 
short length of the chaotic regions. The laminar and chaotic regions of the time series are denoted as Region-I and 
Region-II and are analysed individually. 

We first analyse Region-I of the time series shown in Fig. [T^ & b for linear Granger causality using the procedure 
described in the theory section [ij]. The stress and the VH intensity time series are denoted as X{t) and Y{t), 
respectively. A linear autoregressive model fitting for each of the time series (eqns. [1] and [3]) for increasing model 
order 1 < m < 200 is done and the variance of the prediction errors of ex{t) and ey(t) are calculated. We now 
incorporate information from Y{t) in X{t) and vice versa (eqns. [2]and[4|) and determine the variance of the prediction 
errors of ex\Y{t) ^-nd eY\xit) for 1 < m < 200. The variance of the prediction errors for the mono and bivariate cases 
drops sharply for 1 < to < 30 and saturates for higher values of m. This behaviour is typical of linear regressive 
models and suggests that for the experimental time series analysed here m w 30. Fig. [2^ & b show the error ratios 
var(ej!c|-^(t))/var(ex(i)) and var(ey|jf (t))/var(eY(i)) plotted against the regression model order m. Over the range of 
TO analysed, the error ratios are < 1 implying that X{t) and Y(t) are better predicted when information from Y{t) 
and X{t) are incorporated in the regression equations. The error ratios are at a minimum for m = 33. To understand 
the origin of the minimum in the ratio of the variances, we have evaluated the Akaike Information Criteria (AIC) 
[2^ for various model orders (inset to Fig. [2^). The optimal model order estimated from the AIC yields a regression 
model order m « 34. This suggests that there is maximal driving between the two time series (Fig. [T]) for this value 
of TO and leads to a minimum in the ratio of the variances (Fig. [2]). For this value of m the time series is predicted 
using the various regression coefficients that were computed from eqns. [TK^lEl andU) Figure [3] shows the predicted 
(continuous lines) and the raw (dashed lines) time series of the stress (Region-I of [T^) overlayed on top of each other. 
The predicted time series compares reasonably well with the raw data when information from the VH intensity time 
series is used (Fig. [3}d) than without it (Fig. [3^). Similarly, Figure |4] shows the predicted (continuous lines) and 
the raw (dashed lines) time series of the VH intensity (Region-I of [T|d) overlayed on top of each other. Once again, 
there is reasonable overlap between the predicted and raw time series when information from the stress time series is 
incorporated (Fig. |3|d). This suggests that there exists a feedback mechanism between the stress and the VH intensity 
in the laminar region of the time series. We have carried out the linear causality analysis for the chaotic region of 
the time series also. We find that the variance of the prediction errors do not saturate even for large m (m > 400) 
implying that a linear regression model is inappropriate to model the chaotic regions of the time series. 

In order to determine the influence of nonlinear correlations in the time series on causality, we have carried out 
Extended Granger Causality tests [15]. We describe below the results of this analysis. Following the procedure 
described in the theory section, we first unfold the attractor of the joint dynamics and the attractors for the individual 
time series in a higher dimensional phase space using the delay-embedding technique. The time delay of the two time 
series in the laminar region (Region-I) was found to be r = 7 and was estimated from the first minimum of the 
mutual information function. The embedding dimension was calculated using the method of false nearest neighbors 
and was found to be m = 3. The above two quantities were estimated using the TISEAN software ^2§\. For t — 7, 
the EGCI could not be estimated due to a lack of neighbors, however, t — 1 was found to be sufficient to obtain 
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TABLE I: Shows the various recurrence quantification measures for the laminar and chaotic CRPs 



Recurrence Measure 


Laminar Region 


Chaotic Region 


%Recurrence 


2.2% 


1.7% 


%Determinism 


73.5% 


51.2% 


Maxline 


10 


21 


Entropy 


1.5 


1.4 


Trend 


0.2 


0.6 


%Laminar 


95.1% 


94.5% 


Trap Time 


3.6 


3.9 



causal relations (30| . The error ratios va.T{eY\x{t))/'^^^{^Y{t)) f^nd var(ex|y(t))/var(ex(0) were computed for twenty 
different neighborhoods on the attractor and were then averaged to compute the EGCI defined earlier. The size of the 
neighborhood e was reduced and the above procedure was repeated. Figure [5] shows the EGCI plotted as a function 
of the box size e. We find Ay^x ~ 0.25 a measure of the influence of VH intensity on the stress to be much larger 
than the influence of stress on the VH intensity given by Ax->y ~ 0.02. We are also unable to go to very small e 
due to a lack of neighbors within the neighborhood. We are unable to carry out this test in the chaotic region of the 
time series due to the following reason: an underlying assumption in the above procedure is that are sufficient data 
points within a local neighborhood to perform a linear regression analysis [l5l | . For a high delay-embedding dimension 
m = 7 as in our case, an extremely long time series is required to carry out the Extended Granger Causality test. 

We will now turn our attention CRPs. CRPs were generated and quantified using the software developed by Zilbut 
and Webber [3l| . The stress and the VH intensity time series were normalized to have a zero mean and unit variance. 
The delay and the embedding dimension for the laminar and chaotic regions of the two time series were given in the 
previous section and these were used in the generation of CRPs. The radius of the neighborhood e (fixed radius, 
Euclidean Norm) for the laminar and chaotic regions were 0.38 and 0.61 respectively. These were chosen based on the 
criteria that the %Recurrence should be kept small since large e might cover the entire attractor and all points will 
be found recurrent. Figures [6] & [7| show the CRPs of the VH intensity and stress in the laminar and chaotic regions 
respectively. Inspite of the low %Recurrence « 3% (Table-1) implying a small e, diagonal line structures can be seen. 
These imply that in the laminar and chaotic regions the attractors corresponding to the two time series inhabit the 
same regions of phase space. Not surprisingly, the %Determinism is larger in the laminar compared to the chaotic 
region. The various measures defined above for quantifying CRPs are shown in Table-1. 

V. CONCLUSIONS 

To summarize, we have shown conclusively that there exists a Granger feedback mechanism between the stress and 
the VH intensity time series. The extended Granger causality test, apt for nonlinear dynamical systems, shows that 
VH intensity time series sensitive to orientational order fluctuations has a stronger influence on stress compared to 
the reverse. Most importantly, the extended Granger causality test puts the passive advection approximation used 
in the theoretical model for studying the time evolution of the orientational order parameter on a firmer footing. 
Using CRPs we have also shown that the phase space dynamics of the bivariate time series are similar in the laminar 
and the chaotic regions. 

VI. ACKNOWLEDGEMENTS 

GR's work was supported in part by a grant from DRDO, India and UGC SAP-DSA Phase IV. He and AKS are 
Honorary Faculty Members of Jawaharlal Nehru Centre for Advance Scientific Research, Bangalore, India. 



[1] M. E. Gates and S. J. Candau, J. Phys.: Condens. Matter 2, 6869 (1990). 

[2] N. A. Spenley, M. E. Gates, and T. C. B. McLeish, Phys. Rev. Lett. 71, 939 (1993). 

[3] R. Bandyopadhyay, G. Basappa, and A. K. Sood, Phys. Rev. Lett. 84, 2022 (2000). 

[4] R. Bandyopadhyay and A. K. Sood, Europhys. Lett. 56, 447 (2001); P. Fisher, Rheol. Acta. 39, 234 (2000). 



6 



[5] M. R. Lopez-Gonzalez, W. M. Holmes, P. T. Callaghan, and P. J. Photinos, Phys. Rev. Lett. 93, 268302 (2004); 

W. M. Holmes, M. R. Lopez-Gonzalez, and P.T. Callaghan, Europhys. Lett. 64, 274 (2003); J. -B. Salmon, S. Man- 

neville, and A. Colin, Phys. Rev. E. 68, 051504 (2003). 
[6] A. S. Wunenburger, A. Colin, J. Leng, A. Arneodo, and D. Roux, Phys. Rev. Lett. 86, 1374 (2001); J. -B. Salmon, A. 

Colin, and D. Roux, Phys. Rev. E. 66, 031505 (2002); L. Courbin, P. Panizza, and J. -B. Salmon, Phys. Rev. Lett. 92, 

018305 (2004). 

[7] D. Lootens, H. Van Damme, and P. Hebraud, Phys. Rev. Lett. 90, 178301 (2003). 

[8] R. Ganapathy and A. K. Sood, Phys. Rev. Lett. 96, 108301 (2006). 

[9] A. Aradian and M. E. Gates, Europhys. Lett. 70, 397 (2005). 
[10] S. M. Fielding and P. D. Olmsted, Phys. Rev. Lett. 92, 084502 (2004). 
[11] S. M. Fielding and P. D. Olmsted, Phys. Rev. Lett. 96, 104502 (2006). 

[12] B. Chakrabarti, M. Das, C. Dasgupta, S. Ramaswamy, and A. K. Sood, Phys. Rev. Lett. 92, 055501 (2004); M. Das, B. 

Chakrabarti, C. Dasgupta, S. Ramaswamy, and A. K. Sood, Phys. Rev. E 71, 021707 (2005). 
[13] G. Rienacker, M. Kroger, and S. Hess, Phys. Rev. E 66, 040702(R) (2002). 
[14] C. W. J. Granger, Econometrica 37, 424 (1969). 

[15] Y. Chen, G. Rangarajan, J. Feng, and M. Ding, Phys. Lett. A 324, 26 (2004). 
[16] J. J. Ting, Physica A 324, 285 (2003). 

[17] P. Tass, M. G. Rosenblum, J. Weule, J. Kurths, A. Pikovsky, J. Volkmann, A. Schnitzler, and H. -J. Freund, Phys. Rev. 

Lett. 81, 3291 (1998); M. Le van Quyen, C. Adam, M. Baulac, J. Martinerie, and F. J. Varela, Brain Res. 792, 24 (1998); 

E. Rodriguez, N. George, J. -P. Lachaux, J. Martinerie, B. Renault, and F. J. Varela, Nature (London) 397, 430 (1999). 
[18] S. J. Schiff, P. So, T. Chang, R. E. Burke, and T. Sauer, Phys. Rev. E 54, 6708 (1996); J. Arnhold, P. Grassberger, K. 

Lehnertz, and C. E. Elger, Physica D 134, 419 (1999); R. Quian Quiroga, J. Arnhold, and P. Grassberger, Phys. Rev. E 

61, 5142 (2000). 

[19] M. Casdagli, in: M. Casdagli, S. Eubank (Eds.), A dynamical systems approach to input-output systems, in Nonlinear 

Modeling and Forecasting, Proceedings, vol. XII, Addison- Wesley, 1992, pg. 265. 
[20] S. Boccaletti, D. L. Valladares, L. M. Pecora, H. P. Gelfert, and T. CarroU, Phys. Rev. E 65, 035204(R) (2002). 
[21] J. -P. Eckmann, S. O. Kamphorst, and D. Ruelle, Europhys. Lett. 5, 973 (1987). 
[22] C. L. Webber Jr and J. P. Zilbut, J. Appl. Physiol. 76, 965 (1994). 
[23] J. P. Zilbut and C. L. Webber Jr, Phys. Lett. A 171, 199 (1992). 
[24] J. P. Zilbut, A. Giulani, and C. L. Webber Jr., Phys. Lett. A 246, 122 (1998). 
[25] N. Marwan and J. Kurths, Phys. Lett. A 302, 299 (2002). 

[26] N. Marwan, M. W. Trauth, M. Vuille, and J. Kurths, Climate Dynamics 21, 317 (2003). 

[27] N. Marwan, N. Wessel, U. Meyerfeldt, A. Schirdewan, and J. Kurths, Phys. Rev. E 66, 026702 (2002). 

[28] H. Akaike, IEEE Transactions on Automatic Control 19, 716 (1974). 

[29] The TISEAN pacakage is written by R. Hegger, H. Kantz and T. Schreiber and may be publicly downloaded from 

|http:// www .mpipks-dresden.mpg.de^tisean 
[30] D. Ruelle and F. Takens, Comm. Math. Phys. 20, 167 (1971). 
[31] http: / /homepages. luc.edu~cwebber / ^ (Chuck Webber) 



7 




2500 3000 3S00 4000 4500 

Time (s) 



FIG. 1: Bivariate time series, (a) Stress time series (b) Depolarised Intensity (VH) time series obtained at a shear rate of 




m (model order) m (model ordflr) 

FIG. 2: (a) Ratio of the variance of errors for Y driving X and (b) Ratio of the variance of errors of X driving Y plotted against 
the regression model order m. Here X stand for the stress and Y for the VH intensity. The inset to (a) shows the Akaike 
Information Criteria (AIC) plotted against the regression model order. 
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FIG. 3: (a) Predicted time series of stress (continuous lines) obtained without using information from VH intensity overlayed 
on raw time series (dashed lines), (b) Predicted time series of stress (continuous lines) obtained with using Vlf intensity 
information overlayed on raw time series (dashed lines). Regression model order m = 33. 




FIG. 4: (a) Predicted time series of VH intensity (continuous lines) obtained without using information from stress overlayed 
on raw time series (dashed lines), (b) Predicted time series of VH intensity (continuous lines) obtained with using stress 
information overlayed on raw time series (dashed lines). Regression model order m = 33. 
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FIG. 5: EGCI plotted against the box size e. Y driving X (stars). X driving Y (circles). X stands for the stress and Y for the 
VH intensity. 
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FIG. 6: CRP of stress and VH intensity in the laminar region. 




FIG. 7: CRP of stress and VH intensity in the chaotic region. 



